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We investigate a dynamical Coulomb blockade effect and its role in the enhancement of current- 
current correlations in a three-terminal device with a multilevel splitter, as well as with two quantum 
dots. Spectral decomposition analysis shows that in the Y-terminal system with a two level ideal 
(-y-T ' splitter, charge fluctuations at a level with a lowest outgoing tunneling rate are responsible for a 

r— ^ I super-Poissonian shot noise and positive cross-correlations. Interestingly, for larger source-drain 

f^ , voltages, electrons are transferred as independent particles, when three levels participate in trans- 

^o ' port, and double occupancy is allowed. We can explain compensation of the current correlations 

' as the interplay between different bunching and antibunching processes by performing a spectral 

(-») , decomposition of the correlation functions for partial currents flowing through various levels. In 

the system with two quantum dots acting as a splitter, a long range feedback effect of fluctuating 
potentials leads to the dynamical Coulomb blockade and an enhancement of shot noise. 






o 



X 



PACS numbers: 73.23.-b,72.70.+m,73.23.Hk,73.63.Kv 



o 

^ . I. INTRODUCTION 

For a long time, the interest in current shot noise in nanoscale systems was mostly focused on its reduction from 
y " the Schottky's value Sp = 2e/, where / is the average current, and e is the charge of an electron. The shot noise value 
^_i ■ Sji is reduced because of the Pauli exclusion principle of scattered particles, and for the ballistic point contacts, it 
J^ I can be reduced to zeroi^i^i^ In metallic diffusive wires, the reduction reaches 1/3 x 2elr- and for chaotic cavities, in 
^ . the classic limit, 1/4 x 2eJ— i^, whereas for sequential tunneling through a quantum dot (QD), a maximal reduction is 
' ' 1/2 X 2e/— i^ii^ (for an extensive review see |ll|). 

^ , Electronic correlations can lead to an increase of shot noise above the Schottky's value (so called super-Poissonian 

Q ■ shot noise). lannoccone at ali^ showed that shot noise in a resonant-tunneling diode, biased in the negative differential 

Q [ resistance regions of the I-V characteristic, is enhanced because of the raise of the well's potential energy, which 

causes more states to be available for successive tunneling events from the cathode. Kuznetsov et al.— interpreted 

a transition of shot noise from the sub-Poissonian to super-Poissonian regime in the quantum well as a result of a 

►^ . change of the shape of the density of states, in which a parallel magnetic field leads to multiple voltage ranges of 

^H ' negative differential resistance. Our studies of a ferromagnetic single electron transistor showed^^ an increase of Su, 

00 . much above Sp at the pinch-off voltage; that is, when the current begins to flow, and strong back-scattering leads 

■^ ' to an enhancement in shot noise (see also [lal)- Asymmetry in conducting channels for electrons with opposite spins 

^^ ', can activate spin fluctuations, which results in the super-Poissonian current shot noise.— >^ We call the effect the 

■^ • dynamical Coulomb blockade because the electron that cannot leave for some time the QD blocks the channel for an 

^^ [ electron with the opposite spin. Experimental efforts have just recently been undertake n^^i^^ , in order to search for 

OO , super-Poissonian shot noise in magnetic tunnel-j unctions. The dynamical Coulomb blockade effect can also be seen 

in a nonmagnetic system of two capacitively coupled quantum dots connected in parallel to external electrodes J^ In 

this case, fluctuations of charge polarization can result in an increase of current shot noise. Recent measurements 

were performed on such devices by McClure at al.^'^ and Zhang et al?2i for current noise auto-correlations and cross 

correlations, and confirmed in part the theoretical predictions. Sukhorukov et al.^^ and Thielmann et a\^ suggested 



?H . . . . . 

j^ ■ that inelastic spin-flip cotunneling processes can also lead to super-Poissonian shot noise, which can be observable for 
bias voltages around the corresponding energy for spin-flip excitations. 

In a multi-terminal geometry, one can study the statistics of scattered particles in cross correlation functionsi^ It 
is well known as the Hanbury Brown and Twiss experiment (HBT)^^-, involving two incoming particle streams and 
two detectors showing bunching for photons. The HBT experiments for electrons were performed on a quantum point 
contact^i^ and with fractional quantum Hall effect (FQHE) edge statesSS, and they showed an antibunching effect 
because of the Pauli principle. The cross-correlation in outgoing channels is FQHE effect negative in such case a^^i'^° 
(see also 3l| for cross correlations in a Y-terminal geometry) . Texier and Biittiker showed^ that inelastic scattering 



can lead to positive current-current correlations in a multiterminal system with FQHE edge states, whereas correlations 
remain always negative for quasielastic scattering. In their model, an additional electrode was introduced to keep 
the current equal to zero, and caused voltage fluctuations. Recently, Oberholzer et ali^ verified experimentally 
these predictions in a device, in which interactions between current carrying states and fluctuating voltage were 
controlled by an external gate voltage. Using general arguments for scattering of quantum particles2i, Burkard et 



ali2^ showed that an entangled singlet electron pair gives rise to an enhancement of the noise power (bunching), 
whereas the triplet pair leads to a suppression of noise (antibunching) . Positive cross-correlations were predicted by 
Cottet et al.— i^i^ in sequential tunneling through a three-terminal QD system. The dynamical Coulomb blockade 
effeclii^ii^iii between electrons transferred through a multilevel QD can lead to the super-Poissonian effect for current- 
current auto-correlation functions, as well as to positive cross-correlation functionsj24i2S-36 xhey also analyzed time 
evolution of tunneling events and presented the bunching effect in the system. Gustavsson et al.— ^^ recently took 
time-resolved measurements of electron transport through a multilevel QD, using a nearby quantum point contact 
as a charge detector. They were able to detect bunching of electrons, leading to the super-Poissonian shot noise. 
A phenomenological approach to current cross correlations was presented by Wu and Yip2^. They used a Langevin 
formalism in circuit modeling taking into account voltage and current fluctuations. This general method was applied 
for calculation of current cross-correlation functions in a Y-terminal shaped system. In a similar manner, Rychkov 
and Biittikcr*'' showed that current cross correlations are always positive in a macroscopic classical Y-tcrminal system 
with a fluctuating current in an input branch. 

In this paper, we would like to present a study of a microscopic nature of shot noise in a multitcrminal geometry. 
Our motivation are the recent experimental studies^^'^^ of shot noise on a quantum point contact, acting as a beam 
splitter, and those^- in two capacitatively coupled quantum dots. In the quantum point contact system, the mea- 
surements showed an enhancement of shot noise much above the Schottky value*^ and positive current-current cross 
correlations'*^. Unfortunately, these experiments*^ '^^ do not give an answer to the origin of bunching of electrons. 
In the two quantum dot system, the shot noise results are interpreted^* within a model for a single quantum dot 
with a multilevel structure. The interpretation ignores spatial dynamical charge fluctuations, which occur in tunnel- 
ing through the coupled quantum dots and their influence of shot noise. We expect that the dynamical Coulomb 
blockade effect is responsible for bunching electrons in both the systems, but its origin is different. In a multi-level 
quantum dot, local charge fluctuations are relevant. On the other hand, the potential feedback, caused by spatial 
charge fluctuations, leads to bunching in a multi-dot system. Therefore, our goal is to study quantitatively current 
characteristics and current correlation functions in sequential transport in a three-terminal model with a multi-level 
splitter and with two-quantum dots. 

In the first part of the paper, we will consider the Y-tcrminal system with a two-level splitter showing that super- 
Poissonian shot noise and positive cross-correlations can be expected for large asymmetry in outgoing tunneling rates, 
when a dynamical Coulomb blockade effect can occur. One could naively expect an enhancement of shot noise with 
an increase of a number of levels in the splitter. However, the shot noise is reduced when the third level in the splitter 
becomes to participate in transport. We will show an interplay of auto and cross-correlation partial contributions to 
a total shot noise and its reduction in the Y-terminal system with a three-level splitter. In a multi-terminal system, 
one can change the number of levels participating in transport to different electrodes by applying different voltage 
potentials. We will also present how the current auto and cross-correlation functions can evolve when a voltage 
window is changed in one of the electrode. 

The second part of the paper is devoted to the Y-terminal system with two quantum dots. Our considerations will 
be focused on a feedback effect of potential fluctuations at quantum dots and its role in enhancement of shot noise 
in the current-current auto- and cross-correlation functions. We will show that an origin of electron bunching is in 
this case similar, and it is caused by charge accumulation at one QD. Consequently, this results in strong dynamical 
Coulomb blockade effect. When a second electron participates in transport (for higher bias voltages), antibunching 
processes dominate, and the correlation functions show features typical for sub-Poissonian shot noise. Our analysis 
(in the Appendix) will show that shot noise reduction can be substantial, and the Fano factor can reach its minimal 
value F = 1/(A^ + 1), in a system of N quantum dots connected in series. 

II. Y-TERMINAL SYSTEM WITH A TWO LEVEL SPLITTER 
A. General derivations 

Let us consider a three-terminal system with a two- level QD as a splitter (presented schematically in FiglT]). The 
electronic transport in a sequential regime is governed by the classical master equatio n*'^i'*^i"*'* 
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FIG. 1: Schematic overview of the Y-terminal system with two level splitter, which transfers electrons from the source electrode 
to the two drain electrodes (from the left to the right hand side). 



where po, pi, P2 and p^ denote probabilities of finding an empty system, a system with one electron at the level Ei 
or at the level E2, and a probability of a system with two electrons, respectively. The matrix M is given by 
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where F"™ is the tunneling rate, which changes the occupation of the splitter from the state n to m, F"™ = J2a ^a™'^ 
and F™ = f^nla"^ denotes the total tunneling rate from the a electrode (a = L, Ri, i?2)- 7™" is the net tunneling 
rate through the a tunnel junction, /^^^ = 1/{1 + exp[±(£'„ — iia)/kBT]}, the sign + (or — ) is for a transfer an 
electron from (or to) the electrode, and fia denotes the chemical potential in the a electrode. We assume that the 
energy levels Ei < E2, and when two electrons occupy the splitter, the corresponding energy is E3 = Ei + E2 + U12, 
which includes the electron-electron repulsion energy Ui2. 
The current from the a-th electrode is given by 






(3) 



(The electronic charge in our notation is taken as e.) The probabilities p„ at the stationary state are determined from 
the master equation ((I])-® with the left hand side equal to zero. 

In the sequential regime, it is assumed that tunneling events are independent, corresponding tunneling resistances 
i?a™ ^ Rq — ^/2e^, and the electronic transport is dominated by sequential tunneling processe a'^^''^^ , whereas higher 
order processes (cotunneling) are neglected. Moreover, it is assumed that the resonant current peak for each energy 
level (described by the Brett- Wigner formula) is strongly broadened by temperature; i.e., ^17™' ^ ksT . Comparing 
our results on the shot noise with an experiment, one has to extract a thermal noise contribution, which is always 
present in any conductor (e.g., see [47| for thermal calibration procedure of an experimental setup). 

Fluctuations in the system are studied within the generation-recombination approach for multi-electron chan- 
nels J^i^ The Fourier transform of the correlation function of the quantity X and Y can be expressed asi^ 



Sxy{uj) = 2 / die^"* [{Xmm - {X){Y)] = 4 
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where p„ is the stationary value, X„, and Y^ are the values of X and Y at this state. The conditional probability 
P{n,m;t) to find the system in the state n at time t, if it was in the initial state m at t = 0, satisfies the master 
equation ([TJjiSii^ and its Fourier transform is given by P{n,m;Lu) — [iojl — Af]~^. The elements of the Green's 
function 



G{n, m; to) = [iui — M], 



- Pn/iuj 



(5) 



can be determined directly by matrix inversion. It is useful to perform spectral decomposition, and analyze fluctuations 
corresponding to characteristic eigenfrequencies. Therefore, the Green function ([5]) is represented by its eigenvalues 
Xi, and the corresponding matrix of eigenvectors R as 

G = y^R-'^^R. (6) 

The correlation function between the currents la and la' can be written in the forniiS 

Saa'{uj)=Sa,a'S^''' + S^ao.'{^), (7) 

where 

5f ^'^ ^ 2e(/+ + /- ) = 2e' ^ (rr^™ + r™>„) (8) 

is the high frequency {u; -^ oo) hinit of the shot-noise (the Schottky noise). The frequency-dependent part is expressed 
asia 

s'aa'iu^) = ±2e^J2rr+' ~ C" -^]G™„(c^)[r^ri>„-i - rit'-pn+i] 

+ E[C"+' - r™"-i]G™„(-c.)[rr ^>„-i - rs+^>„+i]} , (9) 

n,m 

where the sign (+) is for the cross-correlation function between the currents in the source and the drain electrode, 
the sign (— ) is for the case with both currents in the drain electrodes or in the source electrode. 

All results for the currents and the correlation functions can be derived analytically using a symbolic mathematical 
program. A problem is presentation of the results, because the formulae can be very complex. Below, we present the 
results for few interesting cases from a physical point of view. 

B. Medium voltage range - Two levels in the voltage window^ 

We consider the case of a moderate source-drain voltage V, for which two energy levels Ei and £'2 lie in the voltage 
window, but the level E3 is beyond the voltage window range and it does not participate in transport. Electrons can 
be transferred only from the left to the right hand side in FiglT] and backflow is ignored. Nonzero elements for the 
total tunnehng rates are: Tf = -fl\ Tf = 7^ , rjo = 7]jO, r|0 = 720 , r^o = 7]^^, rfo = 7^0 . We ignore smearing 
of the Fermi distribution function /^„ and they are taken equal to 1 or 0, respectively, in the expressions for the 
total tunneling rates F^™. In this voltage regime the Schottky term S'f'^'' [Eq.®] is equal to the Poissonian value. 
However, for lower voltages or higher temperatures, when backflow is allowed, S'f '^'* increases and can lead to the super- 
Poissonian shot noise.— li^^i? In this paper we focus only on dynamical processes, which can lead to an enhancement 
of shot noise seen in the term S^^,{ui) [Eq.®]. Therefore, our studies are restricted to the case hj™"- <C fc^T ^ eV. 
In this regime one can find analytical solutions and can easily extract various processes contributing to shot noise. 

Let us begin with a simplest case, when the transfer rates to/from the states Ei and £'2 are equal, i.e. 7°^ = 7^^ = 
7l , 7fl° = 7^" = 7_Ri , 7^*^ = 7^" = 7^2 • The current and the correlation functions at oj = can be expressed as 

. 2e7L7J?.i f.f.. 

27L + iRi + IR2 
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(27L+7fl,i +IR2) 

As one could expect anti-bunching processes dominate in this system. The auto-correlation function Sr^^r^{Q), for 
the currents in the drain electrodes, is always reduced, and the cross-correlation function Sr^^r^ (0) is always negative. 
Now we consider an ideal splitter, i.e. the case when an electron can be transferred only from the level Ei to the 
i?i electrode, but transfer from this level to the R2 electrode is forbidden, and similarly, for transfers from the level 
E2, which can be realized to the R2 electrode only. Thus, the nonzero tunneling rates are: 7^^^ = 7i?i, 7^° = 7/{2 and 



the transfer rates from the source electrode are still kept the same 7^^ 
cross-correlation function are expressed as 



7^^ = 7l. The current, the auto- and the 
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Ax Lo — the current conservation is fulfilled and one can derive all other correlation functions using Eq.(fT3 |) -(fT5 |l . 
e.g. Sr,l{Q) = Sr,rM + Sr,r2{(^). 

Although the tunneling rates 7^?^ and 7^22 ^'^^ different, the currents (|13p in both the drain electrodes are the 
same (as it should be for an ideal splitter). The effect is due to dynamical Coulomb blockade, the process which 
distributes electrons to both the electrodes with the same probability. From Eq. p^ one can determine a condition 
Ihi'^Ri — IR2) > 1r for the super-Poissonian shot noise in the i?i electrode. Positive cross-correlation between the 
currents in the drain electrodes Ri and R2 are for ^l{iri — IR2Y > lRi_lR.2i.lR1 +7-R2)- Both conditions are fulfilled 
for large asymmetry between the outgoing channels, i.e. for tunneling rates 7^?^ ^ 7-R2- 
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FIG. 2: (Color online) Frequency dependence of the Fano factors Fl, Fr^ , Fr^ (black, red, blue curves) - Fig. a, and the current 
cross-correlation functions Srj^r2, Slri, Slr2 (black, red, blue curves) - Fig.b, for a medium voltage window in the Y-terminal 
system with an ideal two level splitter [for the transfer rates 71, — 5, 7^^ = 1, 7^2 ~ 0.5, for which the inverse of the relaxation 
time is 1/r — 0.74 and 10.75 (vertical dashed lines)]. The electron charge is taken as e = —1 in all our presentations. 

In order to have an insight into dynamical processes leading to super-Poissonian shot noise and positive cross 
correlations we perform a frequency analysis of the auto- and cross-correlation functions. Figl2] presents the results 
for a large asymmetry between the tunneling rates in the outgoing channels. The zero-frequency of the Fano factors 
Frj(O), Fl{0) > 1 and the cross-correlations SR^ji2{0), SlRi{,0), Slr2{Q) are positive. It is clearly seen that all 
correlation functions have two contributions corresponding to two dynamical processes with characteristic frequencies 
Ai = 0.74 and A2 = 10.75. For most of the correlation functions the low frequency process leads to an enhancement 
of the auto- and the cross-correlation function, whereas the high frequency process gives a negative contribution. An 
exception is Sl,R2{'^)i for which the lower frequency component is negative and the high frequency contribution is 
positive. It means that various dynamical processes can occur in a system with electron-electron interactions, which 
can lead to bunching or anti-bunching. The bunching process is generally for lower frequency fluctuations, but in 
some cases higher frequency fluctuations can lead to bunching as well. 



C. Three levels in the voltage window 

Now we consider a large voltage regime, for which three energy levels Ei , E2 and that one E^ = Ei + E2 + U12 for 
two electrons in QD, are in the voltage window and all of them participate in transport. The results are presented 
for an ideal splitter, i.e. when an electron can be only transferred from Ei to the i?i electrode, from E2 to the R2 
electrode and transfer from the L electrode to all three levels is with the same tunneling rate. For this case the 
nonzero tunneling rates are: 7^'^ — 7^^ = 7^^, 7^° — ^^ — ^^^ and ^^ = 7^^ = 7^ = 7^ = 7^. The current and 
the correlation functions are expressed by 



Sr^,rA^) = 2e|/flJ 



1- 



Iri = : , (16) 

7L + 7fli 
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Sr,M^) = 0. (18) 



The current Ir^ can be derived from ^TE\i exchanging 7^^ and jr^ , and II = Iri + IR2 ■ Eq- (0 suggests that Ir-^ is 
independent of the tunneling rate jr^ and of the current Ir^ in the second drain electrode. Indeed both the currents 
are independent, which one can see from Ea. (|18p for the cross-correlation function Sr-^_r2{uj) — 0. These formulae 
(fT6|l -(|18 |l are much simpler than those (fT3|) -(fT5 |) for the two levels in the voltage window (when double occupancy of 
the splitter is forbidden) . From comparison of these two cases and one see a role of strong correlations in bunching of 
electrons and enhancement of the auto- and the cross-correlation functions. 

Naively one could expect that shot noise should increase for the splitter with an increase of a number of levels, 
due to an increase of fluctuations through new tunneling channels. A spectral decomposition analysis helps us to 
understand a relation between fluctuations and their contribution to the correlations functions. Using ^ for the 
Green function we can express 



SRi,R^iuj) = 2e|/_Rj 
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where the eigenfrequencies are: Ai = 7l + 7_Ri, A2 = 7l + 7^2; ^'^^ A3 = 2jl + iRi + iR^- The coefBcients are 
s_Ri,_Ri(Ai) = 27^7^^, and sr,^^r^{\2) — sr-^^r-^{X3) = 0. We separate tunneling processes through the level Ei and 

£^3 in the current Ir-^ (t) — Ij^ (t) + Ij^ (t). The coefficients are written in the form 

SR,,R, (A.) = s^^X i^-^) + ^'RuB. (^-) + 4',%i M ' (20) 

where s^™^ (A,y) are the coefficients corresponding to the current-current correlation functions through the lev- 
els n,m = 1,3. Derivations are not complex, and we get s)^j^jj^(Ai) = 1 Rtl\^l l / {l R2 + 7l)^, s^rIpi,S'^2) = 

-iRilR^lllYCiRi +1l){1R2 +7l)], S^W(A3) = iR.lR^llilRi + IR2 +'^lL)|[{lR^ +1l){1R2 + 1l)% SflW(Al) = 

lRilR2ll/{lR2 +1L?, Si?f,^i?i(Ai) = IR1IIKIR2 +lLf- Moreover, the coefficients s^"fl^(A2) = -SrI]rS^2) = 

^R. R (''^2) and s)j ^ (A3) = —Sj^ ^ (A3) = s\^ ^ (A3) for the eigenfrequencies A2 and A3, respectively. This anal- 
ysis shows that for A2 and A3 the inter-level cross-correlation functions have opposite sign to the auto-correlation 
parts and they compensate each other. The nonzero contribution to the shot noise ([TT]) is from the components 
^R R (Ai)i Sr r (Ai) and s]^ ^ (Ai), only. A similar analysis of the cross-correlation function Sr-^^r2{uj) shows total 
compensation of various partial correlation functions for the currents through all three levels. 

D. Different potentials in drain electrodes 

In the Y-terminal system one can apply different voltages to the drain electrodes and get two different voltage 
windows for electron transport. Figl3] presents a situation with a large voltage window for the Ri electrode and a 
medium voltage window for the second drain electrode. Thus, an electron can tunnel from QD to the Ri electrode 
from E3 and Ei. Tunneling to the R2 electrode can be only from E3 (only when two electrons occupy QD); tunneling 
from the single electron state E2 is forbidden because it lies below /i^j. The tunneling rates are assumed to be the 
same as for an ideal splitter, for which the nonzero tunneling rates are: 7°^ = 7^^ = jj^ = 7^ — jl, Jr — Jj^ — iRi , 



/«L 




FIG. 3: Schematic presentation of tunneling for the Y-terminal system with a three-level splitter and with different chemical 
potentials jia in the Ri and R2 electrode. The tunneling rates, which participate in transport are presented. 



7^ = 7^ = 77J2 • The currents and the cross-correlation functions are expressed by 
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The current Ir^ and the auto-correlation function Sr-^,r-^{0) [Eq. pTjl and 123]) ] are the same as for the case of three 
levels in the voltage window [Eq. p^ - P?]) ] considered in the previous section. The cross-correlation function Sr^^^r^ (0) 
is positive for 7^^ > 7^. Performing spectral decomposition ^ one can see contribution of relaxation processes with 
characteristic eigenfrequencies: Ai ~ jl + iRi , ^2 = 1l + IR2 j ^3 = 27^ + ^r^ + ^r^ . The frequency-dependent 
correlation functions are 
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One can see that fluctuations corresponding to the eigen-frequencies Ai and A2 only contribute to the current corre- 
lation functions. A fluctuation process with the characteristic frequency A3 is absent in these correlation functions. 



III. Y-TERMINAL SYSTEM WITH A QUANTUM DOT SPLITTER 
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FIG. 4: Schematic presentation of the Y-terminal system with two quantum dots. The net transfer rates between the quantum 
dots and the electrodes are shown. 



Now, we consider the dynamical Coulomb blockade process in a three-terminal system connected by two quantum 
dots (presented schematically in Figd]). The system is similar to two capacitatively coupled quantum dotsii and an 
experimental setup2i (see the second part of the paper 21] with a Y-terminal structure), where the sub and super- 
Poissonian shot noise were studied. Similarly as in the previous section, we consider sequential electronic transport, 
which is governed by the classical master equationi^'^'^ 
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Here, pa, pi and p2 denote probabilities finding an empty system, a system with one electron either at the quantum 
dot (QDl) or (QD2). The applied voltage difference to the left and the right electrodes is assumed to be moderate, and 
therefore, electron transfers are only from left to right and given by the net transfer rates shown in FiglH Moreover, 
it is assumed that only one electron can be at the splitter, either at QDl, or QD2, and P2 +Pi +Po = 1- 
The currents flowing into the right electrodes are given by 



Iri = ejR^pi = 
Ir2 = e7-R2P2 
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where w = 7^712 -I- 7l7_R2 + 7127^2 + iRilRa- For 712 = 7^^ the current Ir^ = Ir^ and the system of two QD acts 
as an ideal splitter (for any value of jr^ ^^'^ 1l)- 

The current-current correlation functions can be written asi^ 



where the frequency dependent parts are 

52i(^) = 4e2 7iRe[Goi(^)]7LPo, 
S'k.R.i^) = 4e^ 7_R,,Re[Gio(w)]7i?,iPi , 
SR^R^i"^) ^ 46^ 7fl,2Re[G2o(w)]7fl.2P2 
for the auto-correlation functions and 

SlrAuj) = 2e2{7iRe[Goo(w)]7fliPi + lRi'Re[Gii{uj)]-fLPo} , 

Slr2{^) = 2e^{7LRe[Goo(w)]7K2P2 + 7-R.2Re[G2i(a;)]7LPo} , 

SriR2{(^) = 2e^{7_RiRe[Gio(u;)]7i?2P2 +7fl2R-e[G2o(w)]7fliPi} • 



(33) 

(34) 
(35) 
(36) 



(37) 
(38) 
(39) 



for the cross-correlation functions. From Eqs. (|34)l - ((36|) one can see that the autocorrelation functions Saa are propor- 
tional to the current /^ [Ea. pT|) - ([5^ ]. The factor F^ — Saa/'^^Ia is known as the Fano factor. The cross-correlation 
functions Saa> [Eq. ([37|) - ((39|) ] have two terms proportional to the current la and /„', respectively. One can not extract 
a single factor proportional neither to the current la , nor /„/ , nor \/IaIa' ■ It is not possible to define the Fano factor 
for the cross-correlation fmictions [see also Eg. ([43]) for S 11-^112 presented below]. In the literature ^4,49 ,50 (jjffgj-ent 
normalizations of the cross-correlation noise were used and they were called the Fano factor, but in our opinion they 
have no physical meaning. 

At the zero frequency limit the Fano factors for the auto-correlation functions are given by 

Fl = 1 + 2jIji2Jb.i/w'^ - 27L7fl2 [iMlL + 712 + 7i?.i + 7^2) + 7fli7fl.2] /"^^ : (40) 

^fli = 1 + 27L7fli (7L712 - IrJ/w'^ , (41) 

Fr2 = 1- 27L7i27fl2(7i + 712 + 7i?^i + 7^2)/^^ ■ (42) 

In these formulae (|1(I)) - (P^ we extracted positive and negative terms responsible for the super-Poissonian and the 
sub-Poissonian shot noise. From Ea. (|4T|l one can easily see that Fr-^ is in the super-Poissonian regime, if the transfer 
rate jr., is small, and when one can expect a large charge accumulation at QD2. Fig[S^ presents that the Fano factor 
Fl is also larger than unity at small 7^2 • The factor Fr^ < 1 for any transfer rates. It is clear that dynamical 
Coulomb blockade does not occur in the i?2 channel and shot noise is always the sub-Poissonian type. The reduction 
of Fff2 can be substantial, and it can drop to the minimal value Fr^ = 1/3 for 7^ — 712 = jr^ and 7^^ — > [see 
Eq.(|32)]- It is well known that the Fano factor in sequential transport through a quantum dot may be reduced to 
its minimal value F = l/24i2iiiSiii In the present case, however, the system is with two quantum dots connected in 
series. In the Appendix we show that the shot noise reduction can be even larger for a system with a large number of 
quantum dots. The exact calculations show that in a system of N quantum dots connected in series the Fano factor 
can be reduced to the value F = 1/{N -1-1). 

Using the formulae ([57|) - ([55|) one can derive the cross-correlations functions. Below we present the function 

Sr.rA^) = -2IrJr,{jl + 112 + IR, + 2jR,)/w + 2e^jljl2jR,jRjw^ . (43) 

At cj = the current conservation is fulfilled and we can derive all other cross-correlation functions using Eq.([3D])- 
(1121), e.g. SlrAO) = Sr,r,{0) + Sr^r,{0), SlrAO) = Sr,r^{0) + Sr^r^{0). Fig\Ejp presents the results for the cross- 
correlation functions. All of them show a peak at a small jr^ . Increasing the tunneling rate jr^ the correlation function 
SriR2 decreases and changes its sign. It suggests that bunching processes weaken their intensity and antibunching 
becomes dominating. 

Bunching effects and super-Poissonian noise are caused by dynamical Coulomb blockade, which should be seen in 
charge and potential fluctuations in the system (especially at QD2). If we denote a local potential 0i(t) = (t>o + eriiit) 
depending on fluctuation of the number of electrons ni[t) at the i-th quantum dot, then using Eq.([3]) one can write 
the potential-potential correlation function 

^^.^,(0) = 4^0,(n)Re[G„„(O)](/.,(m)p„ = 4e2Re[G,,(0)]ft . (44) 

n,77i 

FiglSt shows that the potential fluctuations S^^if,^ at QD2 are large in the region of small values of 771:2. In FigO: we 
plotted the correlation function Spp for the polarization p{t) = 02 (i) — 4>i{t) between the charges localized at QDl 
and QD2. The corresponding curve (red) presents very large fluctuations in the same region of the parameters as the 
super-Poissonian shot noise. Note that S^-^^^-^ monotonically increases and it does not show activation of any potential 
fluctuations at small 7^2 ■ 

These results suggest that a charge accumulated at QD2 influences the transport through QDl. In order to get 
more information on dynamics of this feedback process we perform a frequency analysis. One can do a spectral 
decomposition of all studied correlation functions and determine components corresponding to eigen-frequencies of 
the system (see eg.Jlq|). For the case presented in FiglHlone gets Ai = 0.5 and A2 = 2.6. Its inverse Ty = l/Xu is 
a relaxation time for eigen-fluctuations in the system. A role of the corresponding relaxation processes in the range 
of the super-Poissonian noise is presented in Fig[6l The Fano factors Fl and Fr.-^ are higher than unity in the low 
frequency range, whereas Fr-^ is always below unity. It means that for F^ and Fr^ their low frequency components 
lead to super-Poissonian shot noise, whereas high frequency contributions reduce shot noise. For given parameters the 
low frequency part can dominate the high frequency contribution and a measurement should show a super-Poissonian 
value of a zero-frequency power spectrum. Fig[6)D shows frequency-dependent plots of the cross-correlation functions. 
The function Sr-^r^ clearly shows two components: low and high frequency ones. Low frequency components are 
positive and they dominate. It means that the low frequency process leads to the positive cross-correlation function 
SriR2 for outgoing electrons. 
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FIG. 5: (Color online) Plot of the Fano factors F^ (for the L-, Ri and R2 electrodes - black, red and blue curve, respectively) 
- Fig.a; the current cross-correlation functions Sac' {SriR2 - black, Slri - red and Slr2 - blue) - Fig.b; and the potential 
correlation functions S^i^i, S^2</>2 (black and blue) as well as the polarization correlation function Spp (red) - Fig.c, as a 
function of the transfer rate 'yR^ at the zero-frequency oj = 0. The other transfer rates are: jl = 712 ~ 7iJi ~ 1- 

FigHt presents fluctuations of the polarization Spp and the charge fluctuations S^2<p2 ^^ QD2, which occur only in 
the low frequency regime. The plot of the function S^-^^-^ is different and shows that the high frequency component is 
also relevant in the potential fluctuations at QDl. Applying the spectral decomposition procedure we can assign the 
lower eigen- frequency to the polarization fluctuations. It means also that the dynamical Coulomb blockade leads to 
the bunching effect for electrons. 



B. Large voltage window 



In Sec lII CI we showed that if two electron states participate in transport through a multi- level splitter, then the 
dynamical Coulomb blockade effect is reduced, and in a special case the cross-correlation function can be Sr-^r-^ {lu) — 0. 
One can expect also a reduction of shot noise for the quantum dot splitter at higher voltages. In this section we consider 
a situation with a large voltage window when energy of second electron overcomes the Coulomb repulsion energy and 
the electron is introduced on an empty QD. The master equation is 





' P12 ' 




' P12 ' 


d 
di 


P2 

Pi 


= M 


P2 
Pi 




. Po J 




L Po J 



(45) 
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FIG. 6: (Color online) Frequency dependence of the Fano factors Fl, Fr-^ , Fr^ (black, red, blue curves) - Fig. a; the current cross- 
correlation functions Sr-^r^, Slri , Slr^ (black, red, blue curves) - Fig.b; and the potential correlation functions 5'0i</)i , 5^2021 
and the polarization correlation function Spp (black, blue and red curves) - Fig.c. The transfer rates are: 7_l = 712 = 7i?i = 1, 
7ii2 = 0.1, for which the eigen-frequency is A = 0.5 and 2.6 (vertical dashed lines). 



where 



M = 



7-Ri - lB.2 


IL 








iRi 


-IL- IR2 


712 


u 


IR2 





-712 - 7fli 


IL 





1R2 


iRi 


-IL 



(46) 



and pi2 is the probability to find simultaneously two electrons at QDl and QD2. The currents flowing into the right 
electrodes are 



'i?i 



e7L7fl.i {IR1IR2 + 7I2 + 7L712 + ILIR2) 
w 

J _ ejLjl2jR2 [iL + iRi + IR2 ) 



(47) 
(48) 



where w = jr^ {-/r^ + jL)ilRi + 7-R.2 + 7l) + 712 blj + 7^2 71- + 7l + 7-Ri (7-R.2 + 7l)]- When the two electron state 
(1, 1) participates in transport the currents increase [compare Eq. (|47|) - (|48)) with Ea. ((3T|l - p2|) ]. but their enhancement 
is different for the different drain electrodes. In previous section we have seen that the system acts as an ideal current 
splitter (with Ir-^ — Ir^) for 712 = 7fli. Now, because single and two electron states are in the voltage window, ideal 
current redistribution is broken (in general, Ir-^ ^ Ir^)- 

The current-current correlation functions are derived analytically in the same way as previously. Because the 
formulae are rather complex, we present plots of the auto- and cross-correlation functions in FiglTl All Fano factors 
are in the sub-Poissonian regime. The factors Fi^ and Fr^ show strong reduction, which can be even below the value 
F = 1/2 for a single quantum dot system. The cross-correlation function Sr-^r^ is always negative. These results 



12 



0.8 



T — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 

t a) 



Ri 



^^ 



0.6 hV"' 



0.4 



0.4 1-"^^ 



§ 0.2 h 



32.. 



I I I I I I I I I I I I I I I I I I I 



^- b) 



-"LR^ 



■~^'^LR| 2 



J I I I I I I I I I I I I I I I I I L 





-0.2 

0.5 1 1.5 2 

FIG. 7: (Color online) Plot of the Fano factors Fa (for the L-, Ri and R2 electrodes - black, red and blue curve, respectively) 
- Fig. a, and the current cross-correlation functions Sa^i (ShiAj - black, Slri - red and Slr2 - blue) - Fig.b, as a function of 
the transfer rate 7^2 at w = and in the high voltage regime, when single and two electron states participate in transport. 
The parameters are the same as in Fig. 5. 

are different than those presented in FiglS] for a mediate vohage regime, when the two electron state is outside the 
voltage window. Now the correlation functions do not show bunching effects - antibunching processes dominate in 
shot noise. 



IV. SUMMARY AND CONCLUDING REMARKS 



Our study addresses dynamical aspects in current shot noise in the Y-terminal system with a multilevel splitter. 
We showed that the zero-frequency current-current auto- and cross-correlation functions increase because of the 
dynamical Coulomb blockade effect. Furthermore, the spectral decomposition analysis indicates that the dynamic 
Coulomb blockade effect is caused by charge fluctuations corresponding to an electron accumulated at the level with 
a lowest outgoing tunneling rate. This process exhibits itself a large asymmetry in the tunneling rates and strong 
Coulomb interactions. In the section II. B, we present the case for the Y-terminal system with an ideal two-level 
splitter, when an electron is transferred from a given level to only one of two drain electrodes, and with a further 
assumption of a single electron occupancy of the splitter. An increase of the voltage changes the situation dramatically. 
For a large voltage, when three levels are in the voltage window and double electron occupancy is allowed, electrons 
transferred to the drain electrodes behave like independent particles with the cross-correlation function Sr^m2 {^) — 
(for any frequency to). Having separated the currents into the partial currents, flowing through each energy level, 
and using the spectral analysis, we were able to decompose the correlation functions, and determine a role of each 
individual component for any eigenfrequency. The partial auto-correlation functions compensate the partial cross- 
correlation functions, leading to Sr^^h2{^) = 0, and reduce the auto-correlation function Sii^^R^{uj) [see Eg. ipTj) ]. For 
the situation considered, bunching and antibunching scattering processes compensate each other perfectly. We also 
considered the case with different voltages applied to the drain electrons, with the voltage windows of different sizes 
(with three and one level in the window, respectively) . The current ffuctuations are then compensated only partially. 
For example, the cross-correlations can be positive and negative, depending on the tunneling rates for incoming and 
outgoing channels. 

We suggest to perform an experiment to see these correlation effects on a multilevel QD system. A direct observation 
of the bunching process on a two-terminal multilevel QD was performed by Gustavsson et al^I^ using a quantum 
point contact system acting as a charge sensor. Our proposal is different because the system should be three-terminal 
one with asymmetric tunnel contacts, and current cross correlation measurements should be performed. In the 
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experiment, one has to change the voltage window, from the low to the strong voltage regime, in order to change the 
number of levels and the number of electrons participating in transport. 

We also considered the Y-terminal system with two quantum dots, which is the model corresponding to the experi- 
mental setup recently studied by Zhang et al.^^. In the experiments^, one can control all the parameters of the system, 
and verify our theoretical predictions. Our studies show that the dynamical Coulomb blockade effect leads to electron 
bunching, and the effect manifests itself with strong charge fluctuations at the side quantum dot (QD2). In order to 
activate the effect, one has to apply the gate voltage, and reduce the tunneling rate ^r^ from QD2 to the electrode. 
We predict an increase of the Fano factor Fr-^ to the super-Poissonian regime and positive cross-correlations Sr-^r^ 
for 7i^2 — > 0. When the bias voltage increases, we expect that the two electron state will become a participant in 
transport, and shot noise will significantly be reduced to the sub-Poissonian regime. We also predict that in a system 
consisting of N quantum dots con (A^ + 1). 




FIG. 8: An experimental setup for two coupled quantum wires, which can be used for studies of shot noise and electron bunching 
(adapted from [5l|]). By tuning voltages applied to the gates G2, G3 and Ge, one can change tunneling transfer rates for the 
incoming and outgoing channels in the splitter, and test the dynamical Coulomb blockade processes in the current-current cross 
correlations. 

We could not point out a microscopic origin of an enhancement of shot noise in the experiments on a quantum 
point contaclj^i^. We also could not find a source of local or spatial charge (potential) fluctuations, which lead to 
bunching in a such simple experimental setup with a quantum point contact. We hope that dynamical Coulomb 
blockade is again responsible for bunching in this case. Therefore, we propose shot noise measurements in a similar 
system presented in FiglBl In this case, two ballistic channels are coupled with each other, and current can be split 
into two electrodes, similar to the experiments in |4lll42l |. This experimental setup allows control of output channels 
by the gate electrodes G3 and Gg, which can also control charge fluctuations in the system. We speculate that 
dynamical Coulomb blockade leads to bunching in coherent electron transport as well. To our best knowledge, there 
is no research in the coherent transport regime showing electron bunching. While there are related papers^ on shot 
noise in the coherent regime, which take into account electron-electron interactions and the Kondo resonance, but 
their shot noise is in the sub-Poissonian regime. 
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APPENDIX: CURRENT NOISE IN QUANTUM DOTS CONNECTED IN SERIES 

We consider a system of A^-quantum dots connected in series in a high voltage hmit, in which electrons can only 
hop from left to right. The master equation can be written as 



dt 



Pn 



P2 
Pi 
PO 



-"fNR JN-IN 






INR 















"723 7l2 

-712 7L1 
-7L1 



Pn 



P2 

Pi 
Po 



(A.l) 



where we have assumed a single electron in the system, i.e. X^n-P" = 1- Here, we present the analytical results for 
the A^ = 3 system. The current is expressed as 



7Ll7l272373fl, 



7-L1712723 + 7il7l273fl + 7Ll72373ii: + 7l272373fl 

The zero-frequency current shot noise is given by 



St,t, — SfiB — St^r — 2el 



LR 



llilhlh + 7ii7?273fl + iIiiIsiIr + ll2lhllR 

(7LI712723 + 7Ll7l273fl + 7Ll72373fi + 7l272373ii:)^ 



(A.2) 



(A.3) 



It is clear that the Fano factor F — S/2eI < 1 and its lowest value min{F} = 1/4 is reached for 7^1 — 712 = 
723 — "f3R- Generalization for any N can be straightforward performed by a mathematical induction, showing that 
min{F} = 1/(A^ + 1). 

This result is a generalization of the Fano factor derived for a single QD in a sequential tunneling regime, when 
min{F} = l/2^ia^ 
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